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. . . Abstract: 

C One of the most important tasks in high energy physics is search for the ex- 

otic states, such as glueball, hybrid and multi-quark states. The transitions 
tp(ns) — ► ip(ms) + 7T7T and T(ns) — ► T(ras) + inr attract great attentions be- 
cause they may reveal characteristics of hybrids. In this work, we analyze those 
transition modes in terms of the theoretical framework established by Yan and 
Kuang. It is interesting to notice that the intermediate states between the two 
gluon-emissions are hybrids, therefore by fitting the data, we are able to deter- 
mine the mass spectra of hybrids. The ground hybrid states are predicted as 
4.23 GeV (for charmonium) and 10.79 GeV (for bottonium) which do not cor- 
respond to any states measured in recent experiments, thus it may imply that 
r S^ ■ very possibly, hybrids mix with regular quarkonia to constitute physical states. 

Oh' Comprehensive comparisons of the potentials for hybrids whose parameters are 

^ ■ obtained in this scenario with the lattice results are presented. 

i-C, . PACS numbers: 12.39.Mk, 13.20.Gd 
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I. INTRODUCTION 

In both the quark model and QCD which governs strong interaction, there is no any fundamental principle 
to prohibit existence of exotic hadron states such as glueball, hybrid and multi-quark states. In fact, to 
eventually understand the low energy behavior of QCD, one needs to find out such states. However, the 
recent research indicates that they may mix with the ordinary hadrons especially the quarkonia. Thus they 
evade direct detection so far, even though many new resonances which have peculiar characteristics, have 
continuously been reported by various experimental collaborations. Theorists have proposed them to be pure 
gluonic (glueball), quark-gluon (hybrid), and/or multi-quark (tetraqurk or pentaquark) structures which are 
different from the regular valence quark structure of qq for meson and qqq for baryon. Since the quark model 
and QCD theory advocate their existence, at least do not repel them, one should find them in experiments. 
However, even with many candidates of the exotic states, so far none of them have been confirmed yet. 
Moreover, the possible mixing of such exotic states with the regular mesons or baryons contaminates the 
situation and would make a clear identification difficult, even though not impossible. From the theoretic 
aspect, one may try to help to clean the mist and find an effective way to do the job. 

The transition of heavy quarkonia such as ip(ns) and T(ns) to lower states tp(ms) and T(ms) (m < 
n) with two pions being emitted, provides an ideal laboratory to study the spectra of hybrids. In the 
transitions ^(ns)(T(ns)) — ► ip(ms)(Y(ms)) + tttt (m < n), the momentum transfer is not large and usually 
the perturbative method does not apply. The QCD multipole expansion (QCDME) method suggested by 

U 3,4,51 



Gottfried, Yan and Kuang[ 



5( well solves the light-meson emission problem. In the picture of the 



multipole expansion, two gluons are emitted which are not described as energetic particles, but a chromo filed 
of TM or TE modes, then the two gluons which constitute a color singlet, hadronizc into light hadrons (|. It 
is worth emphasizing again that the two gluons are not free gluons in the sense of the perturbative quantum 
field theory, but a field in the QCD multipole expansion. It is easy to understand that such transition is 
dominated by the El-El mode, while the Ml-Ml mode is suppressed for the heavy quarkonia case. 

Since two gluons are successively emitted, there exists an intermediate state where the quark-antiquark 
pair resides in a color octet. The color octet q — q and a color-octet gluon constitute a color singlet hybrid 
state. Therefore, in the framework, a key point is to determine the spectra of the hybrid states \qq~g > where 
q can be either b or c in our case. Due to lack of enough data to fix the ground state of hybrid mesons, 
Buchmiiller and Tye 7[ assumed that the observed ?/>(4.03) was the ground state of \ccg >. 



Yan and Kuang used this postulate to carry out their estimation on the transition rates 2|, |3j . For the 



intermediate hybrid states they used the phcnomenological potential given by Buchmuller and Tye 



3 to 



calculate the widths of T(2s) — > T(ls)irn, T(3s) — > T(1s)7T7t,T(3s) — > T(2s)7r7r. The theoretical prediction 
on the rate of T(2s) — > T(ls)7r7r and T(3s) — ► T(2s)7r7r is roughly consistent with data J|, whereas that 
for T(3s) — > T(ls)7T7r obviously deviates from data. It is also noted that when they calculated the decay 
widths, they need to invoke a cancellation among large numbers to obtain smaller physical quantities, thus 
the calculations are very sensitive to the model parameters, i.e. a fine-tuning is unavoidable. Recently 
Kuang [3 1 indicates that determining the proper intermediate hybrid states is crucial to predict the rates of 
the decay modes such as T(3s) — * T(ls)7r7r. 

There have been some models for evaluating the hybrid spectra, but there are several free parameters in 
each model and one should determine them by fitting data. This leads to an embarrassing situation that 
one has to determine at least one hybrid state, and then obtain the corresponding parameters in the model. 
Moreover, the recent studies indicate that hybrid may not exist as an independent physical state, but mixes 
with regular quarkonia states, therefore the mass spectra listed on the data table are not the masses of a 
pure hybrid, which are the eigenvalues of the Hamiltonian matrices. Therefore a crucial task is to determine 
the mass spectra of pure hybrids, even though they are not physical eigenstates of the Hamiltonian matrices. 



1 01 ] collaborations, a remark- 



Recently, thanks to the progress of measurements of the Babar [9| and Belle 
able amount of data on the transitions ip( ns ){Y(ns)) ~~ ^ ip(ms)(Y(ms)) + tttt have been accumulated and 
become more accurate. Since the large database is available, one may have a chance to use the data to 
determine the mass spectra of hybrids. 



In this work, we apply the QCD multipole expansion method established by Yan and Kuang 
potential model given by several groups 



11 



12 



2j and the 



13(, to calculate the transition rates of ip(ns)(Y(ns)) — > 
ip(ms)(T (ms)) + irir by keeping the potential model parameters free. Then by the typical method, namely 
minimizing % 2 for the channels which have been well measured, we obtain the corresponding parameters, 
and then we go on predicting a few channels which are not been measured yet, finally with the potential we 
can determine the masses of hybrids, at least the ground state. 

To make sense, we compare the potentials for hybrids whose parameters are obtained in this scenario with 
the results of the lattice calculation. We find that if the parameters in the potential suggested by Allen et 



al. 



13j adopt the values which are obtained in terms of our strategy, the potential satisfactorily coincides 



with the lattice results. 
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Our numerical results indicate that the ground states of pure hybrid \ccg > and \bbg > do not correspond 
to the physical states measured in recent experiments, the concrete numbers may somehow depend on the 
forms of the potential model adopted for the calculations (see the text). This may suggest that the pure 
hybrids do not exist independently, but mix with regular mesons. 

After the introduction we present all the formulation in next section, where we only keep the necessary 
expressions for later calculations, but omitting some details which can be easily found in Yan and Kuang's 
papers. Then we carry out our numerical analysis in term of the \ 2 method. Comprehensive comparisons 
of various potentials with the lattice results are presented. The last section is devoted to conclusion and 
discussion. 




FIG. 1: 



II. FORMULATION 



A. The transition width 



The theoretical framework about the QCD Multiploc Expandsion method is well established in Refs 2, 
3|, |4j, |5| , and all the corresponding formulas are presented in their series of papers. Here we only make a 
brief introduction to the formulas for evaluating the widths which we are going to employ in this work. 
In Refs. 2j, |3| the transition rate of a vector quarkonium into another vector quarkonium with a two-pion 
emission can be written as 



i 



where \C\ | 2 is a constant to be determined and it comes from the hadronization of gluons into pions, G is the 



phase space factor, f^fif^ i F is the overlapping integration over the concerned hadronic wave functions, 



their concrete forms were given in 3[ as 

,, Pj , Pf _ ^ jRAr)r p -Rh(rydrjR* Kl (r'y P 'Ri(r'y 2 dr' 

Jni,li,n F> l F Mt — Em ' 

where ni,np are the principal quantum numbers of initial and final states, Ii,If are the angular momenta 
of the initial and final states. I is the angular momentum of the color-octet qq in the intermediate state, 
Pi,Pf are the indices related to the multipole radiation, for the El radiation Pj,Pp=l and I = 1. Ri,Rf 
and Rki are the radial wave functions of the initial and final states, Mi is the mass of initial quarkonium 
and Eki is the energy eigenvalue of the intermediate hybrid state. 

B. The x 2 method 

The standard method adopted in analyzing data and extracting useful information is minimizing the x 2 
and in our work, we hope to obtain the model parameters. When calculating % 2 , we would involve as many 
as possible experimental measurements to make the fitted parameters more reasonable. Here we adopt the 



form of x defined in 



1J as 



where i represents the i-th channel, Wf 1 is the theoretical prediction on the width of channel i, W^ xp is the 
corresponding experimentally measured value, AW^ p is the experimental error. 

Wf 1 will be calculated in terms of the potential models with several free parameters which are described 
in the following subsections, thus Wf h is a function of the parameters. By minimizing % 2 , we would expect 
to determine the model parameters. Some details of our strategy will be depicted in subsection E. 

C. The phenomenological potential for the initial and final quarkonia 

In this work, we adopt two different potentials for the initial and final heavy quarkonia and the intermediate 
hybrid states. 

The Cornell potential 



1 51 ] is the most popular potential form to study heavy quarkonia. The potential 



reads as 



V{r)=--+br, (4) 
r 



usually in the literature many authors prefer to use a s instead of n and it has a relation k = ( r ) ; an d 
a s (r) can be treated as a constant for the bb and cc quarkonia. 



G 

The modifed Cornell potential: It may be more reasonable to choose a modified Cornell potential which 
includes a spin-related term [la] , and the potential takes the form 



where the spin-related term V s is, 



V(r) = — + br + V s (r) + V , 



(5) 



with 



5 a (r) = (^Ye 

\/7T 



3„ — cr r 



and Vo is the zero-point energy, ( in Rcf. 



16l | it was set to be zero), here we do not priori-assume it to be zero, 



but fix it by fitting the spectra of heavy quarkonia. 



D. The potential for hybrids 



The intermediate state as discussed above is a hybrid state \qq~g > and we need to obtain the spectra 
and wave-functions of the ground state and corresponding radially excited states. Yan and Kuang used the 
phcnomcnological potential given by Buchmullcr and Tyc 7] to evaluate the mass of the ground state of 
hybrid, instead, in our work, we take some effective potential models which are based on the color-flux-tube 
model. 

Generally hybrids are labelled by the right-handed^^) and left-handcd(n rra ) transverse phonon modes 



and a characteristic quantity A as 



N =Y1 m ( n m + n m)i 



11 



12, 



13 



17 



18 



19] 



All the details about the definitions and notations can be easily found in literature 
Various groups suggested different potential forms for the interaction between the quark and antiquark in 
the hybrid state. We label them as Model 1, 2 and 3 respectively. 
In this work, we employ three potentials which are: 
Model 1 was suggested by Isgur and Paton 



11| as 



V(r) 



h br + -(1 - e" 

r r 



-/6 1/2 r 



Vo- 



(6) 



Model 2: Swanson and Szczepaniak 



12 1 think that the Coulomb term in model 1 is not compatible with 



the lattice results, so that they suggested an alternative effective potential as 

V( r ) = br + -(1- e- fbU ' r ). (7) 
r 

To get a better fit to data, we add the zero-point energy Vq into Eq.0, 

V(r)=br + -(l-e- fbl/2r ) + V . (8) 
r 

Model 3: In model 1, the Coulomb piece is not proper, because the quark and antiquark in the hybrid 
reside in a color-octet instead of a singlet (the meson case) , the short-distance behavior should be repulsive 
(it is determined by the sign of the expectation value of the Casimir operator in octet). Thus Allen etal 



suggested the third model 



131 ] and the corresponding potential form is 



v ( r ) =o+ V(&r) 2 + 2tt& + V . (9) 



Because in these forms the authors do not consider the spin-related term (which we name as V s .), we can 
modify the potential by adding a spin- related term V s , then the potential becomes: 

V(r) = V l + V s . (10) 

By this modification, one can investigate the spin-splitting effects. Generally, V s should have the same form 
as that in ((5|). 

E. Our strategy 

The strategy of this work is that we will determine the concerned parameters in the potential (Eqs.©, 
©, and Eq. (fTU|) ) by fitting the data of heavy quarkonia transitions. 

To obtain the concerned parameters in the potentials (Eqs.©, ©, and (flT)|) ) which specify the hybrids 
sates, we use the method of minimizing \ 2 defined in §Ji§. Concretely, in Eq.Q, W\ h is a function of the 



11[, we set / = 1, therefore x 2 is also a function of those 



parameters k, f, b, Vq and |Ci| 2 , and following Rcf. 
parameters. Minimizing x 2 , one can fix the values of the corresponding parameters. Still for simplifying our 
complicated numerical computations, we choose a special method, namely, we first pre-set a group of the 
parameters, and we calculate the hybrid spectra and wave-functions by solving the Schrodinger equation, 
then we determine |Ci| 2 in Eq.([T]) in terms of the well measured rate of ij)(2S) — ► J/ipirir. With this |Ci| 2 



as a pre-determined value or say, a function of other parameters, we minimize x 2 to fix the values of the rest 
of parameters k, 6, Vo- 

With all the parameters being fixed, we can determine the mass spectra of the hybrids which serve as 
the intermediate states in the transitions of ip(ns)(T(ns)) — > -0(ms)(T(ms)) + tttt. It is noted that the 
spectra determined in this scheme are not really the masses of physical states, unless the hybrids do not mix 
with regular quarkonia. In other words, we would determine a diagonal element of the mass Hamiltonian 
matrix, whose diagonalization would mix the hybrid and quarkonium and then determine the eigenvalues and 
cigen-functions corresponding to the physical masses and physical states which arc measured in experiments. 

III. NUMERICAL RESULTS 

To determine the model parameters in the potential, we need to fit the spectra of ^(ns), r] c (ls), r] c (2s) 
and T(ns) and in this work, we only concern the ground states and radially excited states of cc, bb and ccg, 
bbg systems. 

A. Without the spin-related term V a 

The potentials for quarkonia (Eq.(g])) and hybrid (Eq.© (model 1), Eq.© (model 2) and Eq.© (model 
3)) do not include the spin-related term. In this work, we adopt the Cornell potential to calculate the spectra 
and wavefunctions of the regular heavy quarkonia. The concerned parameters in the Cornell potential have 
been given in literature as for the cc mesons, k = 0.52, b = 0.18GcV 2 , m c = 1.84 GeV, whereas for the bb 



mesons, k = 0.48, b = 0.18GeV 2 , m b = 5.17 GeV 



15J . It is also noted that to meet the measured spectra 



of charmonia and bottonia a zero-point energy Vq is needed. 

The potential for the hybrid takes three possible forms which are shown in Eq.©, ||5J) and ([9]). We keep 
the values m c = 1.84 GeV, rn& = 5.17 GeV which are obtained by fitting the spectra of regular quarkonia 
\bb > (T(ns)) and |cc > (tp(ns)) with the potential (j4|), but need to gain the values of the relevant parameters 
K, b and Vq etc. by minimizing % 2 for the decays ip(ns)(T(ns)) — ► ip(ms)(Y(ms)) + mr. According to the 
measured value for T(ip(2S) — ► J/tpTnr): 

T tot (^(2S)) = 337 ± 13keV 

B{i){2S) -► J/-07T+7T-) = (31.8 ±0.6)% 

B{4>{2S) -» J/i/'ttV) = (16.46 ±0.35)% 
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we express C\ as a function of the potential parameters which exist in the three potentials ( Eqs.©, (JSJ) or 
<j9j) ) and will be determined. It is noted that C\ is a factor related to the hadronization of gluons into two 
pions. so should be universal for both ip and T decays. The parameters in the potentials arc also universal 
for the bb and cc cases except the masses are different. 

Then for T(T(nS) — > T(ms) + tttt) (m < n), wc calculate W\ h in terms of the three potential forms. The 
corresponding experimental values and errors are W i and AW[ xp given in the references which are shown 
in Table [J 



TABLE I: transition rate of T(nS) — * T(ms) ± nn, (in unit of keV) 



decay mode 


Model 1 


Model 2 


Model 3 


Experiment data 


T(2S) -> T(15')7r7r 


9.36 


9.28 


8.69 


12.0 ± 1.8 


T(3S) -» T(15)7T7r 


1.81 


1.67 


1.85 


1.72 ±0.35 


T(3S) -> T(25)7T7r 


0.86 


0.76 


0.86 


1.26 ±0.40 


T(4S) -> T(15)7T7T 


3.87 


3.43 


4.14 


3.7 ±0.6 ±0.7 [10] 


T(4S) -> T(25 , )7T7r 


1.83 


0.2 


1.44 


2.7 ±0.8 [9] 



By minimizing \ 2 ( cc i-©), wc finally get the potential parameters k, b and Vq and the resultant % 2 =4.42 
for model 1, 13.69 for model 2 and 7.26 for model 3 . Then we obtain | d | 2 = 100.39 x 10~ 6 for model 1, 
259.24 x 10~ 6 for model 2, and 121.78 x 10 -6 for model 3, the other parameters are listed in the following 
table (Tabid}. 

TABLE II: potential parameters for hybrid 







6(GeV 2 ) 


Vb(GeV) 


Model 1 


0.43 


0.19 


-0.85 


Model 2 




0.15 


-0.43 


Model 3 


0.59 


0.19 


-0.85 



With these potential parameters, we solve the Schrodingcr equation to obtain the masses of ground hybrid 
states of \ccg > and \bbg > (Tabic W\. It is noted that the resultant spectra depend on the potential forms. 
We will discuss this problem in the last section. 

TABLE III: the mass of hybrids (in units of GeV) 





Model 1 


Model 2 


Model 3 


\ccg > 


4.099 


4.549 


4.226 


\bbg> 


10.560 


11.137 


10.789 
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We also make a prediction on the rates which have not been measured yet (Table IV). 

TABLE IV: prediction(in units of KeV) 



decay mode 


Model 1 


Model 2 


Model 3 


T(4S) -» T(35)7T7r 

1p(3S) -> ^(2S>7T 

ip(3S) -> iP(1S)tyty 


0.60 
14.96 
589.91 


0.57 
14.45 
72.34 


0.61 
14.83 
424.22 



It is noted that values of T(tp(3S) — ► i/j(1S)wtt) predicted by models 1, 2 and 3 are quite apart, while 
r(T(4S) -> T(3S")7T7r) and r(V'(3S") -> ?A(2S*)7r7r) predicted by all the three models are close. 

B. Comparison with the lattice results 



To make sense, it would be helpful to compare the results obtained in our phenomenological work with the 
lattice results which are supposed to include both perturbative and non-perturbative QCD effects. Below 
we show comprehensive com par isons of our potentials with the lattice results. 

3 



Following Refs 



12 



17 



18 



20j . the potentials shown in Fig. 2 are specially scaled by V E +(2ro) which is 
the potential for £+ (N=0) at 2r = 5 GcV -1 (for the vertical axis of Fig.2.). 

In the three graphs of Fig. 2, we present comparisons of the three potentials (models, 1,2 and 3) with the 
parameters fixed in last subsections with the lattice results. In the graphs, the dots are the lattice values 
20 1 . 

It is emphasized that we obtain the potential by minimizing x 2 °f t ne data on ip(ns)(T(ns)) — > 
ip(ms)(T (ms)) + tttt, but do not fit the lattice values. Then our results, especially the third potential 
coincides with the lattice results extremely well. It may indicate that the physics description adopted in 
this scenario is reasonable. It is also noted that by model 1, the short-distance behavior of the potential 
is attractive and obviously distinct from the lattice results. This discrepancy was discussed above that the 
quark-antiquark system in hybrid should be a color-octet and short-distance interaction should be repulsive. 
The second potential (model 2) have the same trend as the lattice results, but have obvious deviations (see 
the graph 2 of Fig. 2). 
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FIG. 2: 



C. With the spin-related terms V s 



3. 



For the regular quarkonia we adopt the non-relativistic potential(NR) Eq.© [16|. Since we add a zero- 
point energy Vq in the potential which can be seen as another free parameter (it is the same for both cc and 
66 quarkonia), we re-fit the spectra of the quarkonia to obtain the corresponding potential parameters in 
Eq.©. We list the resultant values of the parameters in Table fVl In Table IVT1 we present the fitted spectra 
of cc and for a comparison, we also include the results given in Ref. 16( in the table. 



TABLE V: potential parameters for cc 



K 


6(GeV 2 ) 


m(GeV) 


cr(GeV 2 ) 


Vb(GeV) 


0.67 


0.16 


1.78 


1.6 


-0.6 



TABLE VI: Eignvalues for cc in GeV 





JH> 


i>(2S) 


1K3S) 


4(43) 


Vc(lS) 


Vc (25) 


Ref[16j 


3.090 


3.672 


4.072 


4.406 


2.982 


3.630 


this work 


3.097 


3.687 


4.093 


4.433 


2.971 


3.634 



For the 66 quarkonia, the corresponding parameters obtained by fitting data are listed in Table IVHI 
By the parameters we predict = 9.434 GeV, which is consistent with that given by 21 1. 



TABLE VII: potential parameters for bb 
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6(GeV 2 ) 


m(GeV) 


cr(GeV 2 )) 


V (GeV) 


0.53 


0.16 


5.13 


1.7 


-0.60 



Then we turn to the hybrid intermediate states. 

For the hybrids, by the observation made in the previous subsection one can conclude that the third 
potential (model 3) better coincides with the lattice results, therefore, in this subsection when we include 
the spin-related term to discuss spin-splitting case, we only adopt the third potential Eq.([9]). It is reasonable 
to keep the values of m c , mi and a to be the same as that we determined for pure qq quarkonia and we also 
set / = 1. Then following our strategy discussed in previous subsections, we obtain the potential parameters 
which are listed in the following table. 



TABLE VIII: potential parameters for hybrid 





n(ccg) 


n(bbg) 


6(GeV 2 ) 


Vo(GeV) 


the best fitted values 


0.54 


0.40 


0.24 


-0.80 



The fitted values and some predictions are also listed in Tables IIXI and [X] Wc obtain 

|d | 2 = 182.12 x 1(T 6 , 

the mass of hybrids are 4.351GeV, 4.333 GeV for the spin-triplet and spin-singlet cc in the hybrid and 
10.916GeV, 10.913GeV for the spin-triplet and singlet bb respectively. Because of including the spin-related 
term, the "ground states" with the qq (q=b or c) being in different spin structures would be slightly split. 

One can observe that the predicted r(T(45) — > Y(3S)irir) and T(ijj(3S) — ^ ^(2S)tttt) are slightly smaller 
than that predicted in the models without the spin-related term, the future experiments may shed some light 
on it, namely getting better understanding on the mechanisms which one can describe the hybrid structure 
better. 

We also calculate the transition rate of rj c — ► r/ c + tt + it, our result is almost triple that obtained in Ref. 4 1 
and it can be tested by the future experiments. It is noted that since we minimize \ 2 , the decay widths 
that we obtain are different from the central values of the measured quantities. We list the widths wc finally 
obtained in the table IIXI 
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TABLE IX: T transition(in units of keV) 



decay mode 


widths (fit) 


T(25) -> T(1S>7T 


8.73 


T(2S) -> T(15)7T7T 


1.94 


T(35) -> T(25 , )7r 7 r 


0.69 


T(45) -> T(15')7r7r 


4.10 


T(4S) -» T(25)7T7r 


1.88 



TABLE X: prediction(in units of keV) 



decay mode 


widths of predition 


T(4S) -> T(3S>7r 


0.36 


tp(3S) -> i>{2S)ix-K 


8.84 


ip(3S) -> J/ipirn 


12.38 


Tjc(2S) — » ^c7T7r 


335.66 



IV. OUR CONCLUSION AND DISCUSSION 

Search for exotic states which are allowed by the SU(3) quark model and QCD theory is very important for 
our understanding of the basic theory, but so far such states have not been found (or not firmly identified), 
thus it becomes an attractive task in high energy physics. No doubt, direct measurements on such exotic 
states would provide definite information on them, however, it seems that most of the mysterious states mix 
with mesons and baryons which have regular quark structures. Since they are hidden in the mixed states, 
they are not physical states and do not have physical masses, and it makes a clear identification of such exotic 
states very difficult. In other words, they may only serve as a component of physical states. Even though, 
some phenomenological models, such as the color-flux-tube model, the bag model and the potential model 
etc., are believed to properly describe their properties and determine their "masses", in fact, if they mix 
with the regular mesons or baryons, the resultant masses are only the diagonal elements of the Hamiltonian 
matrix. For example, in the potential model, by solving Schrodinger equation, one obtains the eigen-energy 
and wave function, he only gets the element En ~hyb< 4 > \H-hyb\ ( t> >hyb, where the subscript "hyb" denotes 
the quantities corresponding to hybrids. Meanwhile, there is E 2 2 =reg< <P\H reg \<p > reg corresponding to 
the regular quark structure. If the two eigen-states are not far located, they may mix with each other 
and provide an extra matrix element to the hamiltonian matrix, as £12 = E 2 i —hyb< 4>\Hmix\<t> >reg- 
Unfortunately, there is not a reliable way to calculate the mixing matrix element. One may expect to 
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gain definite information about the hybrid states and maybe starting from there he can further study the 
mechanism of the mixing. 

The theoretical framework established by Yan and Kuang confirms that the intermediate states between 
two pion-cmissions in the transition ip(ns)(Y(ns)) — > ip(ms)(Y (ms)) + 7T7t, arc hybrids which contain a 
quark-antiquark pair in color octet, and an extra valence gluon. Based on the color-flux-tubc model, in 
80's of last century Isgur and Paton suggested a potential model for the hybrid, and this greatly simplifies 
the discussion about hybrids and may offer an opportunity to study the regular quarkonium and hybrid in 
a unique framework. After their work, several other groups also proposed modified potentials to make a 
better description on the hybrid states. When Yan and Kuang studied the transitions, there were not many 
data available, i.e. most of the channels were not measured yet. Therefore they assumed that ^(4.03) as the 
ground state of charmed hybrids \ccg > and estimated the transition rates. Thanks to the great achievements 
of the Babar and Belle collaborations, many such modes are measured with appreciable accuracy. Based on 
the experimental data and the theoretical framework established by Yan and Kuang, we minimize the x 
to obtain the model parameters in the potential for hybrid, and with them, we can estimate the masses of 
the ground states of hybrids. The theory of the QCD multi-expansion is based on the assumption that the 
hadronization of the emitted gluons can be factorized from the transition of T(ns)(ijj(ns)) — > T (ms)('ip('ms)) . 
In fact, this factorization may be not complete if the non-perturbative QCD effects are invloved, namely the 
higher twist contribution may somehow violate the factorization. However, as long as the non-perturbative 
QCD effects arc not too strong, this approximation should be acceptable within a certain tolerance range. 
Moreover, in our study, the non-factorization effects are partly involved in the parameter \C\\ 2 of Eq. (1), 
and in our scheme it is also one of the free parameters which are fixed by fitting data. Indeed, it is implicitly 
assumed that \C\\ 2 is universal for all the processes, and it may cause some error. But it is believed that 
since the energy range does not change drastically, the error should controllable. 

In the calculations, we adopt the Cornell potential for the color-singlet qq (q=b or c) system and the 



potentials suggested by Isgur and Paton (model 1) 
Allen et al (model 3) 



111 ], by Swanson and Szczcpaniak (model 2) 



12j and by 



131 ] to deal with the color-octet qq system, we add a spin-related term to the potential 
for hybrid (model 3 only) to investigate possible spin-splitting effects. The numerical results are slightly 
different when this term is introduced. The masses of the ground state hybrids are 4.23 GeV for \ccg > and 
10.79 GeV for \bbg > which are estimated in terms of model 3. When the spin-related term is included, 
the results change to 4.351 GeV, 4.333 GeV for the spin-triplet and spin-singlet cc in the hybrid and 10.916 
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GeV, 10.913 GeV for the spin-triplet and singlet bb respectively. In other two models, the results are slightly 
different. Indeed as aforementioned, a comprehensive comparison of the results with the lattice values, one 
may be convinced that the model 3 may be the best choice at present. All the obtained masses are different 
from the physical states measured in experiments, and it may imply that the hybrids mix with regular 
mesons. 

There are more data in the b-energy range than in charm-energy region. In fact, when we use the same 
method to calculate the transition ip(ns) — > ip(ms) + 7T7t, with n and m being widely apart (say n=4, m=l 
etc.), the theoretical solutions are not stable and uncertainties are relatively large. It indicates that there 
are still some defects in the theory which would be studied in our future works. Moreover, recently Shen 
and Guo 



22j studies the processes in terms of the chiral perturbation theory and considered the final state 
interaction to fit the details of the 7T7T energy and angular distributions. 

The transition of higher excited states of quarkonia into lower ones (including the ground state) without 
flavor change but emitting photon or light mesons is believed to offer rich information on the hadron structure 
and governing dynamics, especially for the heavy quarkonia physics, for example, Brambilla et al. 2s| studied 
the quarkonium radiative decays which are realized via electromagnetic interactions. 

Our studies indicate that the transitions of ip(ns)(T(ns)) — > ■0(ms)(T(ms)) + tt-k may provide valuable 
information about the hybrid structures which have so far not been identified in experiments yet. 

Since we use the method of minimizing x 2 to achieve all the parameters in the potential model for hybrids, 
it certainly brings up some errors. It is a common method for both experimentalists and theorists to analyze 
data and obtain useful information. Definitely, the more data are available, the more accurate the results 
would be. Therefore more data are very necessary, especially the data on the ip families which are one of 
the research fields of the BES III and CLEOc. 
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